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Abstract 

A new method is presented for the simplification of loop integrals in 
one particle irreducible diagrams with large numbers of external lines, 
based on the partial fractioning of products of propagators. Whenever a 
loop diagram in d dimensions has d + 1 or more lines that carry the same 
linear combination of loop momenta, its integral can be reexpressed as 
a linear combination of integrals with no more than d + 1 denominators 
for each such set of lines, of which d are linear in the loop momenta and 
only one quadratic. In multiloop diagrams, the total number of linear de- 
nominators can be reduced further. In integrals with numerator momenta 
there may also be up to d + 1 linear factors in the numerator. 

1 Introduction. 

Tree and loop amplitudes and cross sections with multiple external lines in gauge 
and other field theories are the subject of wide attention [1], [2], [3]. Much of 
this work is based on analyticity properties of the integrals [4] , [5] , and on the 
identification of bases for integrals with fixed number of denominators [6] . 

If the number of loops is kept fixed, one might expect that the diagram 
function of a field theory, in some dimension d, would get progressively more 
complicated as the number of external lines is increased. For one loop, however, 
it has long been known that the number of denominators that are quadratic in 
the loop momentum can be reduced to no more than five for each loop [9] in four 
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dimensions. This paper presents an alternative method for reducing a diagram 
with any number of loops and many external lines to a linear combination of 
integrals of limited complexity. These integrals correspond to vacuum bubble 
diagrams with the same number of loops as the original diagram and with up 
to + 1 denominators per line, d of which are linear in the loop momenta. The 
number of such diagrams increases linearly in L — 1, independent of the number 
of external lines. Such integrals have d + 2 poles per loop, to be compared 
to 2{d + 2) poles in the integral of a standard diagram. The method of Ref. 
[10], [11] can then be applied to derive a further reduction to d denominators, 
at least in the one loop case. The reduction of the number of lines for single- 
particle irreducible diagrams is analogous to the recursive structure of gauge 
theory tree amplitudes, [12] and the underlying graphical identities are in fact 
related [7, 8]. 

The reduction is based on two fairly simple identities, which are presented in 
Section 2. The general formula for the class of what will be called "sufBciently 
dressed" diagrams is also presented there. For such diagrams the total number 
of denominators that depend on loop momenta is k{d+ 1)L, with k a number 
that is defined by the theory and does not grow with d or L. For instance, for 

theory k equals 3. Of these denominators, k x d x L are linear, the others 
quadratic in q. In Section 3 it is shown how to extend the method to theories 
with spinors and vectors. The general result is that there is a basis for the 
integrals such that for each set of lines carrying a specific linear combination 
of loop momenta there is only one quadratic and up to d linear factors in the 
denominator and up to + 1 linear factors in the numerator. It is possible to 
further reduce the number of linear denominators. The method for doing this 
is illustrated by a two-loop example in scalar theory in Section 4. 

2 Identities. 

We consider scalar integrals in d dimensions. To each diagram, we associate a 
vacuum bubble diagram, found by removing all external lines. In this section 
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the method is presented in the case, when the diagram is "sufficiently dressed" . 
A diagram is sufficiently dressed if there are at least d external lines emanating 
from each line in the corresponding bubble diagram. The reduction is performed 
for arbitrary complex masses and complex external momenta. 

In the diagram, consider the set of denominators that corresponds to a given 
linear combination of loop momenta, q. This set of lines can be identified with 
a single line in the corresponding bubble diagram. The starting point of the 
reduction is the following identity, applied to any such set of lines, 

1 1 

n (5 + Pe? +ml= ^J^^^^ {q + Per + ml 

X n - 

X,jj^e ■ ~ + P'j + m'j - pI - ml ' 

where Nq is the number of external lines attached to q, and Pe denotes the sum 
of external momenta carried on line e. This identity is valid as long as none of 
the denominators vanishes. The q,Pe are rf-component momenta, and nie are 
the (possibly complex) masses of the lines. Masses can be zero. The integral is 
regulated by off-shellness. This equation is valid for all Ng > 1. 

Equation (1) is a consequence of the following identity which will be used 
repeatedly below, 

N N 

n ^ = E ^ n ^U=Xe, SeX^ + ae = 0. 

e=l ^ ^ e=l ^ ^ f=l,...,N,f^e •' •' 

(2) 

In this equation, Se,a;, arc arbitrary complex numbers, considered for a 
generic point in the space [sg, ae,x}. The identity (2) is most easily proved by 
considering the integral 

where the contour of integration C encircles zero in a counter-clockwise sense. 
The integral over such contour is equal to the left hand side of Eq . (2). When 
a circle at infinity is added to such a contour, the integration can be closed on 
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all the other poles of the integrand, giving in this way the sum of residues of 
the integrand other than from zero. This is equal to the right hand side of Eq 
.(2). 

In particular, setting Sg = s for all e and taking the limit s — > 0, the identity 
reduces to 

^1^1 1 

e=l e=l /=!,.. .,iV,5te ■' 

Taking fle = (^e +Pe)^ + "Zg gives the identity Eq . (1). 

We next assume that Nq > d. By further partial fractioning the product of 
linear factors in (1) can be further reduced, by using an identity whose proof 
will be given momentarily. For arbitrary Nq> d+1 the identity is 

^ ( n 2n ■ (n- - B ) + I ™2 _ 2 ™2 )^ii,-,j<i(Pe') . 
where the factors ^/i,...,j<j(Pe') are functions of external momenta only, 

In this expression, the d-dimensional vectors Qji,...,j^ are the solutions of the 
following system of d linear equations, corresponding to the "missing" denomi- 
nators in each of the terms on the right-hand side of Eq. (2), 

2g • {Ph - Pi) + pI + ml -pl~ml=0 (7) 
2Q • {Pj^ - Pi) + pI + 'm'ja -Pl-ml=(i ■ 



Thus the argument of Tj^^,,,j^{pe') represents masses and all momenta pe' that 
do not appear in the g-dependent denominators of that term. This sytem is to 
be solved at generic values of momenta and masses, i.e. symbolically. It is worth 
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noting that solutions to Eq. (7) are analogous to the complex momenta that set 
internal lines of the S-matrix on shell in Ref. [12] and discussed diargammatically 
in Refs. [7, 8]. 

The proof of (5) again relies on the identity (2). In order to simplify the 
notation, we rewrite the identity to be proven in the form 

FN{{l},{a},q) = l[-r^—-= (8) 

e=l ■ 9 T "e 



where q e C^, le € C"*, ag G C and 

Te....,e.(W,{a})= n 177^' (9) 

where Q is the solution of the system 

les-Q + aes=0 ,s = l,...,d , (10) 

at a generic point in the parameter space. 

For the proof of (8), single out a particular coordinate system in g-space 
and apply the identity (2) to the component q^, to obtain 

(11) 

where Q'^''{{1}, {a}, Qd-i, qi) is the solution of 1^^ ■ q + = considered as 
an equation for qd- The identity (2) can be applied to the product under the 
sum with qd-i singled out, with the result 

^ 1 

F^{{l},{a},q)=Y^ (12) 

^ E (/ 77+7, l9d=Qa''({'}.W.9d-i,-..,9i)) 

e,_i = l,...,iV,^e, ^'^-^ 9 + «e._i 
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where Q%^''^'^~^ ,s = d,d—l is the solution of the system If ■ q + af = ,/ = 
ed,ed_i, considered as equations for qd,qd-i- Continuing in this way, one may 
arrive at the following identity, valid for any k,l <k < d, 

w 1 

FM{{l},{a},q)=Y,-r^——^---^ (13) 

^ X] , l«r-=Qr'"' ■'°'=+'({0,W,«fc,.-,<Zl), r-=d,..,ft+l) 

^ n (777— r— l9,=Q:<*' - '''=({0,W,9fc-i,---,gi) ,W,...,fe) ' 

where (3^<*' ' ''^''({Z}, {a}, gi), r = d,...,k is the solution of the system 

■ Q + Qes = for Qd,---,<lk- In particlular, for fc = 1 the product under the 
iterated sum is independent of q. It will be denoted by Te^,...,ei({p}, {a})- It 
coincides with the quantity defined in Eq. (9). It is independent of the order of 
the indices e^, ei. Thus, the identity (13) for fc = 1 can be rewritten as 

FN{{l},M,q)= E ^/a,..../i(M'W) 

X E X---X 

{ei,...,ed}={fi,...,fd} 



1 

lei -1 + flei 



where the inner sum extends over permutations of the set {/i, /d}- The 
inner sum is equal to Fd{{lfj^,...,lf^},{af-^,..,af^},q), as can be deduced from 
the identity (13) applied to this particular choice of the momenta and scalars 
Ce, with k — 1, and N — d. This proves the identity (8), and hence Eq. (5) 

Combining the two partial fractionings (1), (5), we see that, for Nq > d, 
the product of Nq > d+1 denominators can be reduced to a sum of terms each 
of which has only d+1 denominators, only one of which is quadratic, 

n {q+pe)^+ml JqTpefT^e ^^^^ 

^ J ^ 29 • {pj^ - pe) +pI + ml -pI- m2^^'i.-'^<^(^^') • 
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We now embed this result in an arbitrary diagram. The integrand of an 
arbitrary IPI Feynman scalar diagram can be written as 

Integrand = g J\ feT^;^^ ' 

Here, i = 1 ... 7 labels the lines of the corresponding bubble diagram and Ei 
labels the set of lines in the diagram associated with line i in the bubble diagram. 
All of the lines in Ei carry the same linear combination of loop momenta, qi, 

L 

Qi = ^ h^cxQa , (17) 

c«=l 

where L is the number of loop momenta. € {0, 1, —1} can be read from the 
bubble diagram in the standard way. 

The partial fractioning formula (15) can be used to derive the following 
algebraic identity for the full integrand of any sufficiently dressed scalar integral 

where Tp^ .i [p^i. ) is obtained through the use of the prescription of Eq. (6) 
with PetPj substituted accordingly. The order of operations is such that right- 
most operations are done first. The argument Pe'. of T's signifies the dependence 
on other momenta from the line i. 

The integrals of both sides of Eq. (18) are well-defined on any representative 
of the homology class in H4l{P'^^ — Ujf?) (D? stand for the standard [4] 
compactifications of the zero sets of the Feynman denominators, and P''^ is 
the complex projective space) that does not intersect the zero loci of the linear 
factors. One such choice corresponds to the standard ie prescription. Therefore 
we may extend the equality of integrands in (18) to an equality of integrals, 



Integral = E + . E , ,^<> («) 

X n 



1 



s=l...d 
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It is enough to choose any representative, C of the homology class and deform it 
in such a way that it does not intersect any of the zero loci of the denominators 
in the reduction formula. Or one can initially choose such a representative. 
There are obviously many such choices. On such representatives the reduction 
formula (18) is valid in each point, so it indeed can be integrated. 
Eq. (18) for the integral can be rewritten as follows 

Integral = V Wt":: [p,,) f (dq)]! \ ^ (20) 

with the single sum over all possible choices of the d-tuples of external lines 
that emanate from each bubble line i, in a sufficiently dressed diagram. The 
rightmost operations are done first. The indices of the T^* j*(Pe^) 

are the same as the indices of the vectors in the second product under the 
integral sign. 

In summary, we see that an arbitrary sufficiently dressed integral can be 
reduced to a linear combination of integrals of simpler type, with one quadratic 
factor and exactly d linear factors for each bubble line. It might be that this 
form of the integral is useful in applications, especially because of the reduced 
number of poles in each integrand. One must keep in mind that there is still 
significant freedom in choosing the contour Q, if one for instance desires to do 
numeric evaluation of this integral, or relate it to physical amplitudes. 

3 Extension to theories with spin. 

In this section it is shown how to extend the reduction method to theories with 
arbitrary spinor or vector numerators. It will be demonstrated that an arbitrary 
integrand can be reduced to a sum of terms such that there is one quadratic 
factor per bubble line, up to d linear factors in the denominator and up to 
d + 1 linear factors in the numerator. Thus, the conclusion that the number of 
elements in this basis for integrals grows linearly with the number of loops and 
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is independent of the number of external lines still holds. 

We start with the observation that for a product of propagators like the left- 
hand side of Eq. (1), the number of quadratic factors in the denominator can 
still be decreased to one through the identity (1) without changing numerator 
momenta at all. It remains to reduce the products of the form 

Ue=l{ie ■ q + ae) 

where g is a loop momentum, Ig are linear combinations of momenta flowing into 
the line and ke are arbitrary momenta, not necessarily related to the Ig- This 
generality of the choice for ke includes open vector indices, so one can choose for 
instance ke,^ = S^,^. The denominator terms tte are expressed through external 
momenta and masses, and the g'-dependent numerator terms be are arbitrary. 
As above, the integer E is one less the number of propagators that correspond 
to the bubble line. Integer r is a free parameter of the formula that accounts 
for the number of linear factors of q in the numerator, ranging from to £J + 1. 
The final case can be realised when all the external lines are attached to the 
bubble line through a vertex that is linear in momenta, such as a triple gauge 
boson coupling. Thus, r = —1, E. 

li r > d then a direct generalization of the scalar formula, Eq. (5), holds 

1= Yl ^ (22) 

with 

where Q is the solution of a system of equations that is analogous to (7) and is 
constructed from the vectors k . Therefore, in this case the basis of intergals 
is the same as in the scalar case. 

In the proof of Eq. (22), the following generalization of the identity (2) will 
be used 

ULlife^ + ^e) ^ 1 , n/^l(//^ + «/) 

ne=l(fea: + ^e) e=l + n/=l,^e + &/) 
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which is valid at a generic point in the space of complex variables fe,ae,ge,be,x, 
and for all integers F < E. The proof involves the same residue argument as 
before and we omit it. 

The proof of (22) now parallels the proof of (5). Consider the quantity 

GfM{1}, {k}, {«}, {b}, q) = nti(fce-g + 6e) ^ ^25) 

where Ig, kg, q G C', Ug, bg € C and F < E — d. Application of the identity (24) 
to the variable qa gives 

GFM{lh {k}, {a}, {b}, q)=Y, (26) 

ed=l • 9 + fled 

, n/=i(fc/-g+^/) , . 

with Q^"* = Q^''{{1}, {a}, qd-i, 3i) the solution for q^ of the equation Ig^ ■q + 
ttg^ = 0. The identity (24) can be applied again to the product inside the sum, 
with qa-i taken as x. The result, analogous to (12), is 

GpMilh {k}, {a}, {b}, q)=Y. , , (27) 

6^ = 1 • 9 + fled 
E 

X 

Ed- 

,Y^akrq + bf) 

n/=l,^e„e,_i(0- •'? + «/) 

with Q^"''"""' = Qr'''''''"'({0;{a},'7<i-2,---,'?i) being the sohition for qd,qd-i of 
the system If-q + at = 0,t ~ ed, ed- 1 • After the application of the basic identity 
(24) s more times the result takes the form 

^ 1 

G^,£({Z}, {fc}, {a}, {6}, 9) = ^ , X • ■ - x (28) 

ed=l 'e<i • 9 + fled 



^ 1 



s + l 



n/=i(fc/-9+M 



^ (l-rE /, , •.\q^=Q'/"""'^-\r=d,...,d-s^- 
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As in the scalar case for s = — 1 can be rewritten as 

GFM{lh{k},{a},{b},q)= E Mf,,...j,{{k},{l},{a},{b}) 

{ei,...,ed}={/i,...,/d} 

with M given by Eq. (23). The sum that multiplies Mf^^,,,j^ is precisely 
Fail e<jj •••5 ^ei J oied) •••) c^ei ) 9)- This provcs the identity (22). 
In case r < d, I can be represented in the form 

1= E ^^/i.-./^ n'^ n\ ^ n + (30) 

/i,...,/<j lL=l(«/a -9 + 01/3) /=B-d+l,...,£-r 

with M the same as above, Eq. (23). This proves the assertion made in the 
beginning of this section. 

Equation (30) can be simplified and recast in the invariant terms, again by 
partial fractioning. We do it here in the case r = d — 1, where for Eq. (30), 

1= E ^/...../.%^^7r^^^ • (31) 

fi,-Jd Us=i{lfs-<l + afJ 

Momenta If^ are linearly independent, so fc^-d+i = J2s'^slfe some ag. 
Therefore, 

1 



(32) 

for some M(<*-i), M^'*). 

It is possible to determine the coefficients M^'^~^\M^'^\ The M^"^^ are given 
by the same formula (23), as can be seen by studying the asymptotics near 
d-tuple denominator poles. It remains to find M^'^~^\ In order to do this, 
consider the line defined by 

Z/i -9 + 0/1=0 (33) 
hd-i ■ 1 + ^fd-i = • 
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Its equation can be written as 



L = {q^ = A^ + B^t\t e C} . (34) 

We multiply Eq. (32) by the factors Z/j -q + af^, ...,lf^_^ ■q + af^_^ and consider 
the result near the line L. It is 

Uf=i,...,E-d+Akf • {A + Bt) + bf) _ ^^^(,_,) ^ Mff...,/, 
n/^A,...,/._,a/-(^ + ^i) + «/) ^^'••'^'^-^ ^lf{A + Bt)+af ■ 

(35) 

It follows from this formula that 

aa^-^) - n/^i,...,i;-d+ifc/ ■ B , . 

Note that B is determined by the d—1 hyperplanes of Eq. (33) up to a constant, 
which cancels in the above fraction. 



4 Further reduction and symmetries for multi- 
loop integrals. 

In this section further partial fractioning in the linear sector of the integrand is 
discussed for the example of two-loop scalar integrals. This allows one to reduce 
the number of linear factors further. Also, a symmetry of the resulting integrals 
is identified. These techniques can be generalized to higher loops. 

In Section 2 it was shown that the scalar integral can be reduced to the 
following sum 

Integral = ^ (pe^)/(pei,Pjj, Sj) , (37) 

where a^i are quadratic functions of masses and momenta that can be read off 
from the previous formulas. In the two-loop case a generic integral I has the 
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form 



liPe. , Pjl ,aei)= j dqidq2 J—^^ 



)2 +mf {q2+Pe2)^ + ml 

n 



(38) 



{qi + 92 +^63)^ + ml ^.^A^i^ . (^^^ - PeJ + dji 



TT TT 



where i G 1,2,3, Ei are d-elemcnt sets of the indices 7'*, and aji arc functions 
of masses and momenta that one can read off from Eq. (10). The change of 
variables 

?i=5l+Pei, (39) 
92 = 92 + Pes , 

will eliminate the dependence of the two quadratic factors on the external mo- 
menta. The prime in q^ will be suppressed in the following. In order to simplify 
the notation, the following change of coordiantes can be made 

r = Pes- Pei - Pe2 , (40) 

Psi = Pji - Pe, , ds- = aj^ - -Pci , Si G , z = 1, 2, 3 . 

In these variables, the integral 7, which will still be denoted by the same letter, 
assumes the form 



n --^ n -— ^ n 



(li-Pai+as,\%. q2-Ps2+as2 (91 + 92 + r) • + 

SiG-C/l S2t-ti2 S3t£/3 

In the general two-loop integral, Eq. (41), there are exactly 2d linear factors 
that have q2 dependence, namely 

^2 = TT — ^ — TT 7 \ • (42) 

p 92 • PS2 + as2 p (91 + 92 + r) • Ps3 + Us^ 

where Ei are d-element sets of indices. The same partial fractioning as in Eq. 
(2) for products of linear denominators can be done on these factors with the 
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result 

(43) 

where we have introduced the subsets E[ and E'2 whose union has d elements, 
as indicated by the first summation over all such subsets. The functions T^, ^, 
are the analogs of the functions T in Eq. (2), 

^^-""^ " n Q'-Ps.+as. n (Q' + gi+r).p«3+a,3 ' ^^^^ 
where Q' is the solution of the following linear system 

Q' ■ + = 0, S2 G (45) 
{Q' + 51 + r) ■ p^3 + as3 = 0, S3 e £^3 . 

The solution of this linear sysytem can be found by inverting the d x d matrix 
whose rows are the d vectors p^^ Pss in Eq. (45), which we assume to be 
linearly independent. It can thus be written 

Q'fi = fi^ + e^uqi,u , (46) 

where /, e are rational functions of all external momenta and are independent 
of qi. The functions T^, ^, in Eq. (44) can thus be written as 

^k,K = n ' (47) 



where F = {E2 - E'^) U (£^3 — £^3) is a rf-element set. We thus have from Eq. 
(43), 

p ^ TT ^ TT 1 TT \ 

|S^|5^^|=d 1^ ■ + '^f' sLe^ «^ ■ + .ai ('^l + 92 + r-) • Pes + a.3 

(48) 

Explicit expressions for v and d are obtained by substituting the solutions (46) 
in Eq. (44). In the resulting form, the 2d g2-dependent denominators P2 in 
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Eq. (43) are transformed to a sum of terms with d linear factors in qi and 
q2, times d linear factors in qi alone. So far, the total number of gj-dependent 
denominators is still the same as in Eq. (43). 

Partial fractioning can now be applied to the 2d remaining denominators 
that depend on qi but not q2, 

= n ,,.'+a ^k,E'. ('9) 
= n n - 

where F is the d-element set in (48). Going through the same steps of partial 
fractioning that led to (43) this quantity can be reduced to a linear combination 
of d linear factors in qi times products of fractions that depend only on external 
momenta, 

where TE'^{p,a) are the familiar T-functions, which in this case are obtained 
from partial fractioning of (49), and the solution of a system of equations, as 
in (45). 

Thus, the original integral /, Eq. (42), can be reduced to linear combinations 
of a new set of integrals with 2d loop momentum-dependent denominators, 
reduced from 3d, 

I{Pei,Pji,aji) = I'{r,ps,,asi) TE[{p,a) , (51) 

\m=d,\E;,\ + \E'^\=d 

where T^j (p, a) is a function of external momenta only and where 

f 1 1 1 

I'ir,ps,,as,) = / dqidq2— ^^-^ 27 i i — 2 (^2) 

J qt + miqi+ [qi + q2 + ry + 

X TT — -^—^ TT — ^ — TT 



^ 91 • Psi + flsi \%, 92 • PS2 + «S2 \%, (91 + 92 + r) • + flss 

Slfc-C/^ S2G-C/2 S3t-C/3 

This looks exactly like the original integrals (41) with the exception that now 
the sets E2 , E!^ in total have d elements (while in the original integral they had 
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2d elements). Quantities m,r,Ps2,as2,Ps3,(is3 are the same as in the original 
integral, while the vectors and the scalars are rational functions of the 
variables of the original integral found by solving linear systems of equations, 
like (45). Note that the intergal /' is still overall UV-convergent, and UV- 
convergent in both qi,q2 sectors. 

Finally, under the change of variables 

q'2 = Q2 + qi+r (53) 
q'i = -Qi 

the function /' transforms as 

where si € E'i,S2 e £'2, S3 & E'^- This identity allows the integral to be reduced 
to a form with no more than [d/2] linear factors with qi + q2 , where [a;] is the 
largest integer not greater than x. Indeed, the integrals I' have the property 
that l-E'21 + I-E3I = d. Therefore, if originally \E'2\ > [d/2], then after such change 
of varibales < [d/2]. 

5 Conclusion. 

This paper describes a new method for the simplification of one-particle irre- 
ducible perturbative diagrams with many external lines. Each set of internal 
lines that carry the same linear combination of loop momentum in any diagram 
is rewritten as a sum of terms with one quadratic and no more that d linear 
denominators in that momentum. In the case of tensor integrals, there can be 
up to d + 1 Hnear factors in the numerator as well. This reduction is achieved 
through the use of two elementary identities, both of which are variants of par- 
tial fractioning, and which are analogous to propagator identities described in 
Refs. [7], [8] for tree diagrams. Arbitrary integrals are linear combinations of a 
basis integrals of this kind, with coefficients that are rational functions of ex- 
ternal momenta and masses, which may take on arbitrary complex values. The 
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analytic continuation of correlation functions is thus simplified. 

The results of this study may offer a new starting point for studies of the 
analytic properties of correlation functions in perturbation theory, related to 
the use of the unitarity method [1],[2]. By writing complex diagrams as sums 
of integrals with reduced numbers of poles, it may be possible to simplify the 
study of both asymptotic behavior and analytic structure. These applications 
will be the subject of future research. 
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